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20.  Abrtract  (Continued) 

triaxial  SQUID  magnetometer  used  to  detect  the  ELF  signal  from  the  Navy  test  transmitter  in 
Wisconsin.  The  required  platform  stabilization  was  demonstrated  in  a  towing  basin  by  motion 
spectrum  measurements  on  a  hydrodynamically  stabilized  buoy  designed  to  be  towed  by  a 
submarine.  Motion  excursions  within  the  ELF  passband  of  30  to  130  Hz  were  about  a  micro¬ 
radian;  hence  motion-generated  noise  must  be  further  reduced  by  about  80  dB.  This  noise  in 
the  ELF  signal  can  be  reduced  by  adaptively  determining  a  vector  approximately  equal  to  the 
earth's  magnetic  field  vector  and  combining  it  vectorially  with  the  I^UID  outputs.  This  adap¬ 
tive  processing  has  been  studied  using  a  computer  simulation  of  the  ^UIDs  and  motion  noise. 
A  prototype  dewar  capable  of  fitting  a  towed  communication  buoy  maintained  a  volume  large 
enough  for  a  triaxial  I^UID  sensor  package  at  the  temperature  of  boiling  liquid  helium  for 
102  days.. 
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DEVELOPMENT  OF  A  SUPERCONDUCTING  ELF  ANTENNA 


INTRODLCTION 

The  Navy’s  proposed  extremely-low-frequency  (ELF)  strategic  communication  system 
is  intended  to  permit  reliable,  secure,  low-data-rate  communications  from  a  single  trans¬ 
mitter  in  the  United  States  to  submarines  at  operational  depth  and  at  ranges  of  several 
thousand  kilometers.  The  skin  depth  of  radio  waves  in  seawater  at  ELF  is  tens  of  meters; 
hence  radio  waves  in  the  band  of  interest  (30  to  130  Hz)  experienced  little  attenuation 
in  penetrating  to  substantial  depths  beneath  the  ocean  surface.  The  required  receiver  sen¬ 
sitivity  in  this  environment  is  approximately  200  dB  below  1  V/m/Hz  (10“^^  T/ZHz), 
a  level  of  performance  that  can  be  achieved  by  cu.nrently  available  E-field  trailing-wire 
antennas. 

To  minimize  operational  constraints  on  submarines  attempting  to  receive  ELF  mes¬ 
sages,  it  is  necessary  to  provide  them  an  omnidirectional  receiving  capability.  The  E-field 
trailing  wire,  which  is  an  elementary  dipole,  has  a  cosine-shaped  azimuthal  antenna  pattern. 
It  is  not  practical  to  design  a  transverse  E-field  sensor  with  adequate  sensitivity  to  com¬ 
plement  the  trailing  wire,  which  means  that  a  trailing-wire  H-field  antenna  is  necessary. 
However  the  best  demonstrated  sensitivity  for  an  H-field  antenna  trailed  at  speeds  of  oper¬ 
ational  interest  is  far  below  the  required  value.  Consequently  we  are  investigating  the 
feasibility  of  using  a  superconducting  quantum  interference  device  (SQUID)  to  provide  both 
the  required  sensitivity  and  the  necessary  omnidirectionality.  As  indicated  in  previous  dis¬ 
cussions  of  this  program  [1,2] ,  a  number  of  engineering  questions  must  be  answered 
before  the  feasibility  of  a  SQUID  for  this  application  can  be  determined. 

One  of  the  most  serious  problems  relates  to  noise  from  rotation  of  the  SQUID  sensor 
in  the  earth’s  magnetic  field  of  5  X  lO"®  T.  Because  the  earth’s  field  is  uniform  on  the 
physical  scale  of  the  typical  SQUID,  a  set  of  orthogonal  sensors  can  be  used  to  sum  the 
vector  components  continuously  and  thus  permit  tne  earth’s  field  to  be  removable  as  a 
dc  component.  One  difficulty  is  providing  enough  dynamic  range  in  the  parts  of  the  sys¬ 
tem  ahead  of  the  vector  summing  stage  to  prevent  the  SQUID  electronics  from  being 
saturated,  and  another  difficulty  is  achieving  orthogonality  among  the  sensors.  If  com¬ 
pletely  free  rotation  of  the  SQUID  in  the  earth’s  field  is  assumed,  the  total  dynamic 
range  requirement  is  about  196  dB.  This  figure  is  far  in  excess  of  the  available  dynamic 
range  of  electronic  components  and,  in  fact,  of  SQUIDs  themselves:  thus  a  certain  degree 
of  motional  stability  must  also  be  provided.  Assuming  that  electronics  can  be  provided 
with  linear  dynamic  range  approximately  equal  to  that  which  is  available  in  SQUIDs,  a 
figure  of  130  to  140  dB,  the  rotation  of  the  sensor  must  be  constrained  to  less  than  10“3 
radian  to  permit  the  necessary  performance. 

These  circumstances  suggest  five  m^or  subproblems  in  demonstrating  the  feasibility 
of  SQUIDs  for  underwater  reception  of  ELF  communications: 


Manuscript  submitted  March  1,  1976. 


1 


DAVIS,  DINGER,  AND  GOLDSTEIN 


•  Providing  adequate  electronic  and  sensor  sensitivity  to  permit  detection 

of  the  ELF  signal  in  noise  of  roughly  T/ZHz; 

•  Providing  130  to  140  dB  linear  dynamic  range; 

•  Constraining  sensor  motion  to  10"^  radian; 

•  Providing  adequate  orthogonality  among  vector  sensors  to  permit  motional 
noise  to  be  removed  by  vectorial  summation; 

•  Developing  a  signal-processing  technique  to  reduce  the  motional  noise  to 
an  acceptable  level  in  real  time. 

In  addition  to  these  five  subproblems,  a  sixth  matter  is  the  need  to  demonstrate  ^at 
a  suitable  cryogenic  environment  can  bo  provided  Tor  a  SQUID,  which  must  be  operatra 
in  a  liquid  helium  bath,  without  unacceptable  logistical  complications,  without  physical 
danger  to  a  submarine  host  vessel,  and  w'ithout  unacceptable  risk  to  submarine  security. 

It  is  the  purpose  in  this  report  to  describe  the  results  of  the  research  that  has  been 
undertaken  to  investigate  these  six  matters. 


MOTION  STABILITY  MEASUREMENTS 

The  noise  environment  of  a  submarine  could  possibly  permit  a  SQUID  to  be  boom- 
mounted  from  a  location  on  the  hull  or  superstructure.  However  submarine-generated 
electromagnetic  noise  may  make  it  necessary  to  tow  the  SQUID  in  a  buoy  as  much  as  100 
meters  away.  In  these  circumstances  the  hydrodynamic  stability  of  the  buoy  will  deter¬ 
mine  the  achievable  limits  on  sensor  motion,  and  such  a  “worst  case”  was  thus  chosen  for 
the  feasibility  testbed. 

An  experimental  buoy  (Fig.  1),  designed  by  the  David  Taylor  Naval  Ship  Rese^ch  and 
Development  Center  to  be  hydrodynamically  stable,  was  subjected  to  extensive  towing 
tests  in  the  David  Taylor  towing  basin.  The  buoy,  which  is  3.7  meters  long  and  0.5  meter 
in  diarr.eter,  was  instrumented  with  a  triaxial  rate  gyroscope  and  a  triaxial  linew  acceler¬ 
ometer  and  was  towed  in  a  650-meter  basin  at  between  2  and  7  knots.  Several  different 
body  trim  angles  were  investigated  at  each  speed,  and  runs  were  repeated  several  times  to 
insure  statistical  significance  of  the  data.  The  top  trace  in  Fig.  2  is  a  recording  of  the  out¬ 
put  from  one  rate  gyroscope  as  an  example.  The  time  integral  of  this  output  is  shown  in 
Fig.  2  (bottom  trace)  in  order  to  display  the  actual  motional  excursions  in  milliradians. 

An  oscillation  approximately  0.4  milliradian  in  amplitude  with  a  period  of  2  seconds  is 
evident.  The  larger  oscillation  with  a  period  of  about  20  seconds  is  an  artifact  of  the  signal 
processing. 

The  data  were  also  analyzed  by  a  fast-Fourier-transform  technique  to  display  the 
spectral  composition  of  the  noise,  and  examples  of  these  noise  spectra  are  shown  3 

for  towing  speeds  of  4  and  7  knots.  The  data  show  low-frequency  motions  of  the  body 
that  increase  in  bandwidth  as  speed  is  increased.  The  4-knot  spectrum  occupies  a  band 
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Fig.  1  —  Towed  communication  buoy  used  for  the  motion  spectrum  measurements 


RUN  31 

PITCH  SCNSOR 
7  KNOTS 


Fig.  2  Real-time  output  of  the  pitch  rate  gryoscope  (top 
trace)  and  the  time  integral  of  the  top  trace,  displaying 
pitch  angular  excursion  as  a  function  of  time  (bottom  trace) 
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Fig.  3  —  Spectre  of  the  pitch-rate  gyro¬ 
scope  output  for  two  towing  speeds 


of  about  0  to  2  Hz,  and  the  7-knot  bandwidth  extends  from  0  to  4  Hz.  There  «  "o 

below  10“®  radian  in  the  0  to  10  Hz  band. 

Tvoical  accelerometer  data,  which  provide  information  at  higher  frequencies  are 
showl^rrt  5  The  data  are  expressed  in  radians/s^  (angular  acceleration)  relative  to 
10-3  riianfs2  The  spectral  peaks  at  35  and  42  Hz  are  due  to  vortex  shedding  from 
Wunt  trailing  edges  of  the  buoy  hydrofoils  and  could  be  reduced  by  proper 
rru  oolr  at  fin  H7  is  due  to  DOwer  frequency  contamination.  Integration  of  these  data 

th/ 30  U,  130  H.  ELF  bandwidth  can  be  con- 

strained  to  less  than  lO"®  radian. 

To  summarize  the  motion  sUbOity  mea«.mmen«, 
are  10-3  ,adlan  at  most,  are  the  largest  in  amphtude  and  hence 

rtHTLoToS'.o°“  Tetba^rLTlO-V^dtrglmtes  a 
siimri  of  10-11  T  which  is  76  dB  above  the  receiver  sensitivity.  The  signal-processing 
Snique  that  will  ^  dt^  1.  capable  of  reducing  the  in-band  motion  noise  by  at 

least  this  amount. 

Further  detaUs  on  the  motion-spectrum  tests  can  be  found  in  Ref.  3. 
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Fig.  4  ~  Integral*  of  the  spectra  shown 
in  Fig.  3,  displaying  pitch  angular  ex¬ 
cursion  on  a  logarithmic  scale  a*  a  func¬ 
tion  of  frequency.  In  each  case  the 
sensor  resolution  curve  displays  the 
computed  sensor  noise  level. 


FREQUENCY  (Hr) 


Fig.  6  —  Spectrum  of  the  x-oriented  ac¬ 
celerometer  output.  The  spectrum  noise 
floor  can  be  observed  between  80  and 
100  Hx. 


CRYOGENIC  ENVIRONMENT 

Present  technology  can  solve  the  engineering  problems  of  storing,  handling,  and  cycling 
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simplest  and  most  reliable  method  of  providing  the  required  cryogenic  environment,  for  a 
SQUID  that  may  have  to  be  towed  in  a  buoy  100  meters  from  its  host  vessel,  would  be  to 
equip  the  buoy  with  a  container  that  can  carry  enough  liquid  helium  to  keep  the  sensor 
cool  for  a  typical  maximum  submarine  deployment  period  of  about  90  days.  Some  of  the 
considerations  that  contributed  to  this  decision  are  described  in  Ref.  4. 

Accordingly  an  investigation  was  made  of  the  feasibility  of  storing  in  a  container 
small  enough  to  fit  into  a  typical  communication  buoy  enough  liquid  helium  to  provide 
90  days  cooling  for  a  SQUID  system  of  appropriate  dimensions  to  satisfy  the  ELF  com¬ 
munication  requirement  discussed.  This  investigation  required  the  design  of  a  horizontal, 
superinsulated  dewar  container  with  a  180-liter  liquid-helium  capacity.  Figure  6  is  a  draw¬ 
ing  of  the  dewar  that  was  designed  by  Superconducting  Technology,  Inc.,  to  meet  these 
requirements.  It  is  a  right  circular  cylinder  1.8  meters  long  and  0.56  meter  in  diameter 
and  will  fit  into  an  AN/BRA-8  communication  buoy.  Figure  7  is  a  photograph  of  the 
dewar  that  was  procured  to  test  the  dewar  design.  This  dewar  successfully  maintained  a 
SQUID  at  liquid  helium  temperatures  for  102  days  and  thus  demonstrated  the  feasibility 
of  providing  the  necessary  cryogenic  environment.  The  dewar  performance  is  further 
discussed  in  Ref.  5. 


SENSITIVITY,  DYNAMIC  RANGE,  AND  ORTHOGONALITY 

Requirements  on  sensitivity  and  linear  dynamic  range  have  been  described.  Physical 
orthogonality  of  the  sensor  elements  is  also  related  to  an  adaptive  process  that  will  be 
described  for  the  removal  of  noise  from  motion  of  the  triaxial  SQUID  sensor  in  the  earth’s 
magnetic  field.  In  particular,  if  the  sensors  can  be  made  orthogonal  to  an  adequate  degree, 
then  certain  error  terms  in  the  adaptive  signal  processing  will  be  small  and  will  greatly 
simplify  the  computations  that  must  be  made  in  the  adaptive  process. 

A  triaxial  SQUID  sensor  was  procured  for  test  purposes  with  the  specifications  listed 
in  Table  1.  The  device  was  designed  and  manufactured  by  SHE,  Inc.,  and  employed 
shielded  point-contact-type  SQUIDs  that  were  coupled  to  pickup  loops  wrapped  on  a  pre¬ 
cisely  machined  quartz  block  about  7  cm  on  a  side.  Figure  8  is  a  photograph  of  the  pick¬ 
up  assembly  and  shielded  SQUIDs.  The  grooves  in  the  quartz  block  accommodate  both 
the  triaxial  pickup  loops  and  a  calibration  coil.  Leads  from  this  assembly  are  coupled 
through  the  fiberglass  mounting  to  the  SQUIDs,  which  are  in  a  doughnut-shaped  super¬ 
conducting  shield  immediately  behind  the  mounting.  Figure  9  is  a  photograph  of  the  pick¬ 
up  loop  assembly  and  shielded  SQUIDs  with  a  protective  can  in  place,  showing  the  mount¬ 
ing  boom,  baffles,  and  radiofrequency  circuitry.  This  entire  assembly  is  placed  into  a  liquid 
helium  bath,  with  the  SQUIDs  and  pickup  coils  at  the  bottom.  The  three  SQUID  axes 
have  independent  cables  to  control  circuits  that  can  be  placed  as  far  as  100  meters  from 
the  dewar. 

The  triaxial  test  SQUID  has  been  tested  in  the  noise  environment  of  the  surface  at¬ 
mosphere,  which  is  approximately  40  dB  noisier  than  its  design  environment,  with  an 
electromagnetic  shield  made  of  superconducting  lead  around  it  to  permit  its  sensitivity  to 
be  checked  at  lower  noise  level,  and  also  on  the  ocean  floor  at  a  depth  of  100  meters. 
Figure  10  is  a  comparison  of  ELF  transmitted  signals  from  the  Navy  test  transmitter  in 
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overall  length  183  meters 
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Fig.  6  ~  Long-hold-time  dewar  designed  to  fit  a 
communication  buoy 


Fig.  7  —  Long-hold-time  dewar  having  the  design  shown  in  Fig.  6  and  used  to  test  the 
feasibility  of  storing  enough  liquid  helium  to  cool  a  SQUID  for  90  days 
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Tiitilf  I  SjHi  lfication!.  of  the  .SQUID  KI.K- Antenna  Test  Device 


Parameter 

(ioal 

Aehieved 

Sensitivity 

-20n  dB  re  1  V'/m/Hz 

-198  dB  re  1  I  'm/Mz 

;  Linear  ilynamie  range 

140  dB 

139.2  dB 

Orthogonality 

10"**  radian 

0.7  X  lO"**  radian 

.3KiB  bandwidth 

30  to  130  Hz 

0  to  130  Hz 
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Ki«.  9  —  Photo(r*ph  of  th»  SQUID  antenna  probe.  The  flux  tranaformer  aaarmbly  shown  in  Fig  8 
w  hidden  by  the  protective  fiberglasa  shield  on  the  left-hand  side  of  the  probe.  (The  scale  is  in 
inches.) 


F*.  10  —  Comparison  of  transmitted  ELF  signals  received  on  the  SQUID 
antenna  and  on  a  conventional  air-core  loop  antenna.  The  asimuthal 
loop  and  azimuthal  SQUID  traces  are  offset  vertically  for  display 
purposes. 
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Wisconsin  as  received  on  the  triaxial  SQUID  and  on  a  conventional  air-core  loop  antenna. 
Two  of  the  three  SQUID  axes  appear  and  are  shown  independently,  because  in  this  sta¬ 
tionary  test  no  motion-noise  processing  was  necessary.  The  azimuthal  signal  is  the  main 
component  of  the  transmitted  energy,  but  at  receiving  locations  near  the  transmitter  a 
radial  component  also  exists.  This  radial  component  is  shown  for  illustrative  purposes 
on  one  SQUID  axis  in  Fig.  10  but  is  unsealed.  It  is  apparent  that  the  signal-to-noise 
ratios  of  the  loop  and  SQUID  azimuthal  channels  are  similar  (the  SQUID  and  loop  chan¬ 
nels  are  offset  on  the  figure  to  permit  convenient  comparison  —  their  actual  signal  levels 
are  the  same),  and  that  the  temporal  behavior  of  the  signals  received  on  them  is  qualita¬ 
tively  similar.  Slight  differences  in  frequency  response  and  signal  processing  between  the 
two  systems  are  responsible  for  the  small  differences  in  the  appearance  of  the  two  signals. 


SIGNAL  PROCESSING  FOR  MOTION-NOISE  REDUCTION 

The  towing  tests  described  earlier  indicate  that  buoy  motion  in  the  30  to  130  Hz  band 
can  be  constrained  to  be  less  than  10“®  radian.  However  additional  motion  compensation 
must  be  employed  to  reduce  the  level  of  motion-induced  noise  still  further.  The  signal  in¬ 
duced  in  a  magnetic -vector  sensor  rotated  through  an  angle  D  in  an  earth’s  magnetic  field 
from  an  initial  orientation  Oq  of  the  sensor  axis  relative  to  is 

H„  =  //g[cos  00  (cos  n  -1)  +  sin  0o  (1) 

For  Hg  =  5  X  10-5  T  and  H  =  lO"®  radian,  in  a  1-Hz  bandwidth  the  resulting  is  76 
dB  above  the  sensitivity  requirement  of  10-^“*  T/JHz.  In  this  section  an  adaptive  signal¬ 
processing  technique  will  be  discussed  which  can  provide  at  least  this  remaining  degree  of 
motion  compensation. 

An  adaptive  technique  is  required  because  tbe  output  of  a  SQUID  sensor  is  indepen¬ 
dent  of  the  absolute  value  of  the  magnetic  field  along  the  sensor  axis.  That  is,  the  SQUID 
sensor  follows  external  magnetic  field  changes  accurately,  with  a  low  noise  level  and  true  dc 
response,  but  the  actual  magnetic  field  value  along  the  axis  of  the  sensor  is  never  known. 
However  this  magnetic  field  value  must  be  determined  to  remove  motion  noise.  The  es¬ 
sence  of  the  technique  used  to  serve  this  purpose  is  that  by  sampling  the  changes  in  mag¬ 
netic  field  that  occur  as  the  SQUID  sensors  undergo  motion  in  the  earth’s  field,  and  know¬ 
ing  the  functional  form  which  these  changes  mu.st  assume,  then  the  absolute  values  of  the 
magnetic  field  components  can  be  determined  by  the  parameter  estimation  methods  that 
have  been  developed  for  adaptive  control.  The  adaptation  technique  that  has  been  devel¬ 
oped  for  this  purpose  is  described  more  fully  in  Ref.  6  and  will  be  summarized  here. 

For  sensors  that  are  perfectly  orthogonal  a  combination  of  the  three  sensor  outputs 
may  be  formed  that  is  independent  of  sensor  orientation.  This  combination  is  simply  a 
vector  sum.  If  is  defined  as  the  sensor  output  due  to  motion  in  the  earth’s  dc  field 
Hg,  the  quantity  I  l^is  invariant  and  equal  to  f/| ,  because  in  the  absence  of  signal 

and  noise  there  is  no  other  source  of  magnetic  field  strength  than  H^.  Thus 

|2  =  2H^  •  +  H2  +  h2  =  Hf  (2) 
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2H,  •  H„  +  h2  -  0 


(3) 


In  addition  to  an  output  generated  by  motion  in  the  uniform  field  of  the  earth,  the  SQUIDs 
also  respond  to  the  ELF  signal  vector  H,  .  so  that  the  total  sensor  output  is 


Ho=H„  +H 


Computing  2H^  '  Hq  +  H§  and  making  use  of  Eq.  (3)  gives 


2He  •  Ho  +  Hg  =  2H, 


H,  +  2H 


m 


H.  + 


(4) 


(5) 


.  1.L  j  I  LI  I  in-l4  r  a  lower  limit,  and  the  maximum  in-band 

'll"  J^om'h  "f  U  W-S  T  Substitution  of  these  values  into  Eq.  (5)  indicates  that  the  second 
thW  ted  o;“he  dt^nd  side  are  neariy  dx  Oiders  of  ma^itude  lower  than  the 
first  term.  Hence  there  remains  an  error  quantity  S  defined  as 


S  =  2Hw  •  Ho  +  Hg  =  2H, 


(6) 


fhe  quantity  S  is  proportional  to  the  signal  .ector  and,  when  properly  e^luated,  is  inde- 
jendent  of  sensor  orientation. 

The  ineiusion  ‘jj' ^""^"ents 
f/h*  LTmodiM  by'^husu  deviations  from  the  ideal  ease;  these 

H  are  relabeled  and  become  the  components  d],  42,  and  .43  of  a  vector  ^ 

that  is  only  approximately  equal  to  H.  .  The  equation  lor  S  takes  the 

form 


S  =  2A  ■  Ho  +  H' 


(7) 


A  is  used. 

The  adaptive  processing  then  proceeds  as  follows:  ^he  SQUID  outputs 

"iSi  r 

put  proportional  to  H,.  This  final  output  is  free  of  motion-related  noise. 

An  error  analysis  shows  that  for  the  processing  to  reduce  motion  "oi*® 
the  deviations  of  the  sensor  axes  from  orthogonality  in  radians  and  the  2 

tivity  errors  of  the  SQUID  outputs  must  each  be  less  than  radian  as  mea- 

where  as  before  S2  is  the  peak  platform  motion  in  radians.  If  D  -  10 

eiirpH  in  the  motion  stability  tests  described,  then  this  quantity  has  the  value  2  X  10  . 

Table  1  shows  that  the  deviations  of  the  sensor  axes  than 

0.7  X  10-4  radian;  and  the  fractional  sensitivity  errors  of  the  SQUID  putputs  are  less  man 

10-6. 
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The  principal  task  of  the  adaptive  processing  is  to  determine  the  values  of  the  com¬ 
ponents  of  A  rapidly  and  accurately.  Two  techniques  have  been  investigated  to  determine 
the  adaptive  vector  A .  The  first  is  a  least-squares  technique  that  determines  A  by  mini¬ 
mizing  the  integrated  value  of  over  a  time  interval  Tq.  The  second  is  a  continuous 
steepest  descent  method  that  updates  the  current  value  of  A  according  to  the  integrated 
value  of  S. 

These  two  adaptive  techniques  have  been  investigated  using  a  computer  simulation 
(Fig.  11)  of  the  SQUID  sensors  and  the  adaptation.  The  input  signals,  labeled  0,  in  Fig. 
11,  are  the  tape-recorded  outputs  of  the  three  rate  gyroscopes  that  were  mounted  on  the 
buoy  during  the  motion  stability  measurements  described  earlier.  The  outputs  of  the  rate 
gyros  are  the  angular  time  rates  of  change  of  the  buoy  orientation  in  three  orthogonal 
directions.  Integration  of  the  rate  gyro  signals  gives  the  angular  excursions  of  the  buoy. 
By  combining  the  buoy  excursions  with  assumed  values  of  the  components  of  the  earth’s 
magnetic  field  in  an  Euler  transformation  matrix,  magnetic  field  signals  are  generated  that 
simulate  the  motion-induced  magnetic  field  variations  that  would  have  been  observed  if 
SQUID  sensors  had  been  mounted  in  the  buoy.  The  effects  of  nonorthogonal  sensors  are 
included  by  a  transformation  from  an  orthogonal  to  a  nonorthogonal  coordinate  system 
whose  deviations  from  orthogonality  are  specified  by  the  angles  These  procedures 
yielded  the  inputs  to  a  computer  simulation  of  the  least  squares  fit  and  the  steepest  de¬ 
scent  processing  scheme.  The  output  of  the  adaptive  processing  S  was  Fourier  analyzed 
by  a  standard  fast-Fourier-transform  (FFT)  program.  The  simulated  SQUID  outputs  were 
also  Fourier  analyzed  for  comparison  with  the  Fourier  transform  of  S. 


Fi?.  11  —  Computer  simulation  of  the  adaptive  proceuing.  The  initial  ettinutes  are 
required  only  for  the  steepest  descent  technique. 


i 
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U  io  pnmnuter  simulation  was  to  establish  that  the  least  squares 

ties  and  thereby  remove  motion  noise  above  30  Hz. 

The  results  of  the  processing  using  the  data  from  one  typical  lowing  tertmnwUl  be 
f  A  The  sneed  was  7  knots  and  the  run  duration  approximately  90  seconds,  al 
ZStnly  .pprLltely  30  Mcond,  of  the  run  could  be  proceeaed  at  one  time  t»c.u» 

of  computer  limitations. 


Least  Squares  Adaptation 

a  *  a.  K-ia  19  t trace  A\  disolays  a  SQUID  output  (Vj)  as  an  example  of 

theln7:tatSre“e::it^ao*ua^« 

H.,  were  held  conattnt  of  one 

E ?hii‘eZ?s^5ronto1  darwtth'’L“ 

feHrl‘‘SHHS55!S%r“ 

SuHSS  = 

value  of  2He  '  H,,  which  is  10-18  f2_ 


Steepest  Descent  Adaptation 

The  steepest  descent  adaptation  technique  also  was  investigated.  It  was  found  that 
The  steepest  aescem  a  p  materially  affecting  S,  but 

'  l^ruTthe^on  Ji^  to  ‘hnh  »“«*  value,  and 

S'^rnraCrh: 

™rrhthi:srrhrs^»:ud^^^ 

angles  was  also  similar. 
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Fig.  12  —  Computer  simulation  results  for  the  least  squares  technique.  Trace  A  is  the  time  variation  of 
the  SQUID  output  Vj .  Traces  B  and  C  display  the  time  variation  of  the  estimated  values  of  Aj  for  a 
value  of  Tq  of  2  and  15  seconds  respectively.  The  exact  value  of  A  ^  is  also  shown  on  the  left-hand  scale 
in  both  traces.  Trace  D  is  the  time  variation  of  the  error  quantity  S. 
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Fig.  13  —  Frequency  ipectrum  (top)  of  the  SQUID 
output  Vi,  g.ven  by  trece  A  in  Fig.  12,  «nd  fre¬ 
quency  ipectri  in  (bottom)  of  the  error  quantity  S, 
given  by  trace  D  in  Fig.  12. 


The  simulation  demonstrates  that  the  least  squares  technique  ® 
rat«  estimates  of  the  adaptive  constants  without  the  need  for  any  initial  “tirnatc  ^  How 
ever  the  method  is  not  a  closed-loop  technique  and  requires  up  to  10  . 

evaluate  the  adaptive  constants,  which  is  a  decided  handicap,  ^he  requu-em  . 

microorocessor  might  be  considered  a  disadvantage  in  some  applications,  but  in  the  sub 
mSrcomlnTcation  system  under  consideration  here  a  microprocessor  can  Pro^bly 
orovided  The  essentially  continuous  determination  of  the  adaptive  constants  offer  y 
Sie  steepest  descent  technique  and  its  closed-loop  nature 

decent  tX^^  superior,  because  it  can  be  implemented  by 

nutinc  machinery  On  the  other  hand  the  requirement  of  an  initial  value  of  each  of  the 

SapUvrrom  Jts,  which  must  apparently  be  rather  close  to  the  true  value,  ts  a  drsadvan- 

tage. 

An  optimum  system  might  consist  of  a  hybrid  combination  of  the  two  techniques 
Processing  could  be  inlttated  by  aeroing  the  SQUID  outputs,  entering  10  second,  rt  .,Q.  .u 
^^lo  a  computer  memory,  and  soiling  for  the  adaptive  constants  usmg  the  least 
sauLes  technique.  This  value  would  then  be  used  to  begin  the  mtegration  of  the  .  oeepest 
Ascent  eouations  which  could  be  by  analog  or  digital  techniques.  Continuous  monito.mg 

of  the  output  S  by  a  threshold  detector  would  allow  ''«j"^ID  XlluSd  be' 
in  the  earth’s  field.  If  the  noise  in  S  became  excessive,  the  SQUID  outputs  could  ne 

zeroed  and  the  determination  begun  anew. 
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CONCLUSIONS,  FUTURE  PLANS,  AND  RECOMMENDATIONS 


It  has  been  shown  that  the  cryogenic  environment  necessary  for  SQU  D  use  in  a  sub¬ 
merged  receiving  application  can  be  provided  simply  and  reliably  by  suitable  packaging  o 
a  liquid-helium  dewar.  A  horizontal  dewar,  of  an  appropriate  size  to  fit  m  a  typical  com¬ 
munication  buoy,  has  been  procured  and  tested  and  has  demonstrated  a  102-day  hold  time. 

Motion-stability  measurements  have  been  made  on  a  hydrodynamically  stabilized  buoy 
in  a  towing  channel  at  the  David  Taylor  Naval  Ship  Research  and  Development  Center  and 
have  shown  that  motion  in  the  C  to  30  Hz  band  can  be  constrained  to  lO'^  radi^.  Mo¬ 
tion  within  the  30  to  130  Hz  ELF  reception  band  was  held  to  less  than  10  radian. 

In  view  of  these  estimated  buoy  performance  limits,  a  triaxial  test  SQUID  was  pro¬ 
cured  with  a  139-dB  linear  dynamic  range,  a  value  that  will  permit  the  SQUID  ^  be  used 
aboard  a  towed  buoy  without  electronic  saturation  due  to  motion  of  the  SQUID  in  the 
earth’s  magnetic  field.  The  SQUID  also  was  configured  with  a  triaxial  pickup  assembly  in 
which  the  three  pickup  loops  were  orthogonal  to  within  0.7  XlO'  radian.  This  high  de¬ 
gree  of  orthogonality  simplifies  the  adaptive  signal  processing  scheme  that  is  necessary  to 
employ  a  SQUID  in  this  application. 

Because  a  SQUID  cannot  sense  the  absolute  value  of  the  earth’s  magnetic  field,  the 
earth’s  field  strength  must  be  determined  adaptively  to  a  degree  of  precision  that  wi  1  per¬ 
mit  it  to  be  removed  from  the  vector  sum  of  the  sensor  outputs  to  below  noise  level.  Two 
adaptive  techniques  were  studied  and  evaluated  with  the  aid  of  simulated  motion  data 
from  the  motion  stability  measurements.  A  hybrid  of  open-loop  (least  squares)  and  closed- 
loop  (steepest  descent)  adaptation  methods  was  found  to  be  the  most  attractive  approach 
to  solving  the  signal-processing  problem. 

The  principal  remaining  effort  in  this  exploratory  development  program  will  be  de¬ 
voted  to  the  following  matters. 

•  Perform  off-line  adaptive  processing  tests  to  evaluate  adaptive  methods 
using  a  PDP  11/40  minicomputer; 

•  Procure  high-speed  analog/digital  electronics  to  interface  the  SQUID 
directly  to  the  computer.  Perform  further  adaptive  processing  tests. 

•  Measure  the  ELF  noise  spectrum  of  a  submarine  in  the  San  Clemente  test  range; 

•  Tow  a  SQUID  system  in  a  buoy  and  record  the  ELF  signal  from  the  Navy 
test  transmitter.  Evaluate  the  adaptive  signal  processing  (off  line). 

These  investigations  will  permit  the  feasibility  of  actual  SQUID  operation  in  the  under 
water  environment  to  be  demonstrated.  They  wUl  also  permit  on-line  adaptation  to  be 
demonstrated  in  principle.  However  engineering  of  an  entire  communication  i^eiver  and 
signal  processor  for  operation  in  an  underwater  environment  will  not  be  undert^en.  fur¬ 
thermore,  the  adaptive  processing  method  will  not  be  refined  to  the  point  at  which  true. 
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•11  u  .iKio  Rnth  of  these  m&ttBfs  £Uf6  left  until  a  further  inves- 
real-time  adaptation  environment  is  completed.  Some  preliminary  obser- 

tigation  of  the  underwater  ELF  conducted  as  the  opportunity  arises  during 

attractive. 

The  following  Item,  are  recommended  for  the  adwmeed  development  program; 

•  Submarine  noise  survey, 

*  •  Boom-mount  development, 

•  Towed  buoy  development, 

•  Real-time  signal  processor  development,  and 

•  Field  verification. 

The  boom-mount  and  develop- 

a  clear  preference  is  e*Pre“ed  by  '®^  °' should  take  2^ears  and  should  cost 

ment  programs.  The  advanced  *^^heduling  delays  and  field  platform  in- 

approximately  $1  mUlion,  exclusive  of  field-test  scheduimg  aeiay 

stallation  costs. 
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